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ABSTRACT 

We present the reddening (E(V — I)) and Extinction maps in y-band (Ay) and /-band 
(Aj) for 48 Optical Gravitational Lensing Experiment II (OGLE-II) Galactic bulge 
(GB) fields, covering a range of —11° < I < 11°, with the total area close to 11 square 
degrees. These measurements are based on two-band photometry of Red Clump Giant 
(RCG) stars in OGLE-II VI maps of GB. We confirm the anomalous value of the 
ratio of total to selective extinction Ryr = Ay/E( V — I) = 1.9 ~ 2.1, depending 
on the line of sight, as measured bv lUdalskl ((2003) . By using the average value of 
Ryi — 1.964 with the standard deviation sdev — 0.085, we measured E(V — I), 
Ay and Aj, and we obtained extinction and reddening maps with a high spatial 
resolution of 26.7" ~ 106.8", depending on the stellar density of each field. We assumed 
that average, reddening corrected colours of red clump giants are the same in every 
field. The maps cover the range 0.42 < E(V - I) < 3.5, 0.83 < Ay < 6.9 and 
0.41 < A] < 3.4 mag respectively. The zero points of these maps are calibrated by 
using V — K colours of 20 RR Lyrae ab variables (RRab) in Baade's window. The 
apparent reddening corrected /-band magnitudes of the RCGs change by +0.4 mag 
while the Galactic coordinate I varies from +5° to —5°, indicating that these stars are 
in the Galactic Bar. The reddening corrected colour of RRab and RCGs in GB are 
consistent with colours of local stars, while in the past these colours were claimed to 
be different. 

Key words: dust, extinction - Galaxy:bulge - Galaxyxenter - stars:horizontal branch 
- stars:variables:othcr 



1 INTRODUCTION 

A study of stellar populations and stellar dynamics in the 
Galactic Bulge (GB) is important for understanding how 
bulges formed, what are their populations, gravitational po- 
tential and structure. 

Several gravitational microlensing survey groups have 
found hundreds of events towards the Galactic center 
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and thousands are 



expected in the upcoming years by MOA I , OGLE-III 2 
and other collaborations. The data from such microlensing 
surveys is useful to study the Galacti c structure by mea- 
suring the microlensing optical depth |^daEk^^alj^^93 
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suited for numero us other scientific projects (see IPaczvriskil 
Il996t lGouldlll996r) . 

However, as is well known, the extinction due to the 
dust is very significant towards the GB. This affects the 
Color Magnitude Diagram (CMD) of the field and makes 
a separation of stellar populations difficult. To correct for 
these effects the measurements of the extinction in these 
fields are crucial. 

ISchleeel. Finkbeiner fc Davisl dl99ct) made all sky ex- 
tinction map using COBE/DIRBE data, which overestimate 
the extinction tow a rds the GB because of bac kground dust 
jDutra et al.l2003h . lSchultheis et alJlll999ll and lDutra et alj 
(2003) constructed if -band extinction maps of the Galactic 
central region with a resolution of 4' by using J and K 
photometry of the upper giant branch stars in DENIS and 
2MASS data respectively. Some determinations of the ex- 
tinction towards Baade's Window have been performed with 
a number o f different tech niques, including those of stellar 
simulation <Ng et alJll996ft. mean magnitudes of red-clump 
stars (iKiraea. Paczvhski fc Stan ek|l997l ), the a bsolute mag- 
nitude of RR Lvrae stars nfAl^ocket^nil998bl) and magni- 
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tude of the K-giants llGould. Popowski fc Terndrudll99sl) . 
The Large-Scale Extinction Map of the Galactic Bulge was 
made by using the m ean colour of all stars in the MA CHO 
Pro ject Photometry |p" opowski. Cook fc Beckei)l2003f) 

IWozniak fc Stanekl il996T> proposed a method to in- 
vestigate the ratio of total to selective extinction based 
on two-band photometry of Red Clump Giants (RCGs). 
The RCGs are the equivalent of the horizontal branch 
stars for a metal-rich population, i.e., relatively low-mass 
core helium burning stars. RCGs in the Galactic bulge oc- 
cupy a distinct reg ion in the colour magnitude diagram 
flStanek et alJ [2000 and references therein) . The intrinsic 
width of the luminosity and colour distribut ion of R CGs 
in the Galactic bulge is small, a bout 0.2 mag dStanek et alJ 
ll997l:IPaczvnski fc Staneklll998h . 

The CMD is used to obtain the quantitative val- 
ues of the offset on the CMD between the different sub- 
fields, caused by differential extinction. They used RCG- 
dominated parts of the CMDs for determining the offsets, 
the clump being seen at fainter magnitudes and redder 
colours in subfields with higher extinction. They applied this 
method to the OGLE-I data and then found the ratio of to- 
tal to selective ext inction Ryi = A y IE{V — I) = 2 .44. This 
is consistent with iNg et alJ Jl996ft . IStanekl ( 1996) applied 
this method to the OGLE-I data to obtain differential ex- 
tinction Ay and reddening E(V — I) in a 40' x 40' region of 
Baade's window, w ith resolution of ~ 30" . The y estimated 
R VI = 2.49 ± 0.02. iPaczvrlski et alJ il999ft and ISumi eta!] 
l2003al applied this method to OGLE-II (14'. 2 x 14'. 2 with 
resolution of 20" x 20") and MO A (16 deg 2 with resolution 
of 3.45' x 3.45') data respectively. They first made a redden- 
ing map for their fields because determining the reddening 
E(V — I) (horizontal shift in the CMD) is easier than Ay 
and Ai (vertical shift in the CMD). Then the extinction 
map was calculated according to the following formulae: 

Ai = RixE(V-I), (1) 
A v = R V iXE(V-I). (2) 

with "standard" values of Rvi = 2.5 and /?/ = Ai/E(V - 
I) = 1.5. 

iPaczvnski fc Stanekl il998Tl and 

IStutz. Popowski fc Gouldl 1^999) found that the mean 
V — I colours of GB RCGs and RR Lyrae, dereddened 
with Stanek's map, are redder than colours of their nearby 
counterparts. iPopow ski (2000) summarized possible expla- 
nations of these discrepancies, and noted that the simplest 
and the most plausible explanation is a non-standard 
interstellar extinction. The discrepancy would v anish if 
Ryi = 2.1 rather than the standard value of 2.5. lUdalskj 
(2003) showed that there is indeed an anomaly in the 
extinction law, with Ryi = 1.9 ~ 2.3, depending on the 
direction of the line of sight. 

In this paper we confirm the anomalous value of Rvi, 
and by using new value, we construct extinction maps for 48 
Galactic Bulge fields obser ved by the Optical G ravitational 
Experiment 3 II (OGLE-II: ludalski et al.ll2000D . 

In §1^1 we describe the data. We measure the reddening 

3 see |h ttp://www.astrouw .*edu.pl/~ogle| or 
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line in §[3] We construct the reddening and extinction maps 
in §2] Discussion and conclusion are given in §0 



2 DATA 

We us e the VI photometr ic maps of standard OGLE tem- 
plate (lUdalski et al]l2002f) . which contain VI photometry 
and astrometry of ~ 30 million stars in the 49 GB fields. 
Positions of these f ields (BUL_SC1 ~ 49) can be seen in 
lUdalski et ai] |2002ft . We do not use BUL.SC44 in this work 
because most of RCGs in this field are close to, or even below 
the IZ-band detection limit of OGLE due to high extinction. 
The photometry is the mean photometry from a few hun- 
dred measurements in the /-band and several measurements 
in V-band collected during the second phase of the OGLE 
experiment between 1997 and 2000. Accuracy of the zero 
points of photometry is about 0.04 mag. A single 2048 x 
8192 pixel frame covers an area of 0.24 x 0.95 deg 2 with 
pixel size of 0.417 arcsec/pi xel. Details of the instrumenta- 
tion s etup can be found in lUdalski. Kubiak fc Szvrna riski 
<1997t) . 



3 THE RATIO OF TOTAL TO SELECTIVE 
EXTINCTION 

To measure the ratio of total to selective extinction, i.e. 
Ri = AijEiy — I), we make use of the position of RC Gs in 
the (I,V - I) CMD, as it was done bvlUdalskil fcOOft. bu t 
contrary to IWozniak fc Stanekl jl996n and IStanekl <199d) . 
who used the (V ,V — I) CMD. The reddening corrected I- 
band magnitude of the RCGs does not vary with colour, 
while the l^-band magni tude is a function of V — I 
iPaczvnski fc St anek||l99 8l). and using it can lead to sys- 
tematic errors dPopowskil 12000ft . We make a preliminary 
assumption that the average, reddening corrected RCG 
colour is constant, {V — I)o.rc = 1.0, which is approxi- 
mat ely the colour of nearby R CGs as measured by Hippar- 
cos dPaczvnski fc St anck 1998). The mean colour (V — /)rc 
and /-band magnitude (/)rc of RCGs follow reddening line 
with a slope Ri and a constant /o,r,c, to be determined for 
every field 

</}rc = /o,rc + Ri((V - /)rc - (V - I)o,rc). (3) 

In this Section, we measure the slope Ri and brightness 
of RCGs Io t Rc in Eq[3] with a low spacial resolution, but 
with a significant number of RCGs. In the next Section, we 
measure the colour of RCGs with high spacial resolution by 
using these Ri and Io,rc- This is because an identification 
of RCG centroids in / is more difficult than in V — I because 
of the vertical structure of Red Giants overlapping RCGs in 
the CMD. 

In this Section, there are 2 steps to find the RCG cen- 
troids in the CMD. At a first step, we divide the field into 
"bins" to take spacial differences of the extinctions into ac- 
count, then we measure the rough mean colour of RCGs in 
each bin. At a second step, we combine bins with similar 
RCG colour into group to enlarge the significance of RCGs, 
then we estimate the RCG centroids by the Gaussian fit. 

As a first step, we divide each field into 16 x 64 = 512 
"bins" which have 128 x 128 pixels each. As these "bins" are 
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relatively large there may be a considerable differential red- 
dening within them, and the RCG may be elongated along 
the reddening line. 

We select an elongated window within which the RCGs 
are located, following the reddening line as given by the Eq. 
J3J| with a width of ±0.5 mag in /. As we do not know the 
correct reddening slope and the correct magnitude of RCGs, 
we adopt a broad range of trial values: 0.6 < Ri < 2.0, and 
14 < Iq^rc < 15. We also adopt a broad range of colours for 
the search of RCGs, selected to be within 1.4 < (V — I) < 
5; the boundaries are adjusted for each field to minimize 
the contamination by blue disk main sequence stars and 
faint red bulge main sequence stars. We measure the average 
(V — /}rc in each bin using 2-a clipping. These values are 
used as the initial values of the RCGs colour in the next 
paragraph. 

In a given bin we measure the average colour and mag- 
nitude of RCGs within the circle with a radius of 0.4 mag in 
the CMD centered at the colour (V — 1}rc on the redden- 
ing line. These new average values: (V — I)rc and {I)nc, 
calculated for all 128 x 128 pixel bins, are used to obtain 
an improved value of Ri and Io,rc- This process is iterated 
until the values of Ri and Iq,rc do not change any more. 
We found that this final value for a given OGLE field was 
independent of the first guess in the ranges 0.6 < Ri < 2.0 
and 14 < Io,rc < 15, and are roughly the same as the more 
precise values measured in the second step. 

Even if we locate the window slightly higher or lower, 
i.e. 7o,rc in Eq.@, than the true RCG centroid in /, we 
would get roughly same colour as the true RCGs colour be- 
cause the RCGs colour are very similar to the colours of red 
giants which are somewhat brighter or fainter than RCGs. 
We use only resultant average colour {V — I)rc of each bin 
in the following analysis, to arrange bins in order of their 
extinction in a given field. So, as long as we fix the value of 
Ri and Io,rc in a given field these colours can be used to 
arrange the bins. 

As a second step we arrange the bins in a given field in 
order of extinction by using (V — /)rc- Then we combined 
these bins into groups from low extinction to high extinction 
until each group is filled by ~ 1000 RCGs. Given a large 
number of RCGs in each group of bins we could find the 
RCGs positions in the CMD for each group independently. 
We measured RCGs centroid of each group in the CMD by 
following three methods: 

[1] In a given group of bins we measure the average 
colour and magnitude of RCGs within the circle with a ra- 
dius of 0.3 mag in the CMD, centered at the initial values of 
colour and magnitude of RCGs. These new average values 
are used for new RCG selection. This process is iterated un- 
til the value {V — 7)rc and (/)rc do not change any more. 
We found these final values do not depend on their initial 
values which can be given by the rough positions of RCGs 
in previous section, as long as their real positions are within 
this initial circle. 

[2] Similar with [1] but with a larger radius of 0.45 mag 
and weighting V — I and I with two dimensional Gaussian 
with o = 0.15 mag. 

[3] We fit the distribution of stars in each group with 
the power law plus Gaussian luminosity function: 
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Figure 1. Sample /-band Luminosity Functions of stars in rela- 
tively low (upper panel: BUL_SC1, grp=3) and high (lower panel: 
BUL_SC5, grp=10) extinction fields. Solid lines are the best fits 
by equation J1J. Dashed and dot-dashed lines represent the Gaus- 
sian and power law components in the best fits. 

MI) = Pol0 pi/ + P 2 exp (- (J ~J/^ C)2 ) , (4) 

where po, pi, pi and ai t nc are free parameters, calculated 
for each group of bins separately. RCGs are selected within 
the circle with radius of 0.4 mag centered at this best fit 
{I)rc and (V — I)nc measured by method [1]. We fitted the 
colour distribution of selected RCGs with another Gaussian: 

4>vi{V -I) =p 3 exp I — ^ '—^ '- — '—], (5) 

where p$ and ovi, rc are free parameters. 

In Fig.^we show the sample /-band Luminosity Func- 
tions of stars in a relatively low (upper panel: BUL_SC1, 
grp=3) and high (lower panel: BUL_SC5, grp=10) extinc- 
tion fields with the best fit by equation J1J. Here groups 
(grp) which have roughly 1,000 RCGs are numbered from 
1 in order of increasing extinction. One can see how signif- 
icantly we can measure the RCGs centroids in these large 
groups. 

We arrange 44 OGLE-II fields into 11 regions of fields 
close to each other, as presented in table U Region s A, B, 
C and D are identical with those in Udalskjl J20031) . Fields 
BUL_SC6, 7, 47, 48, 49 are not analyzed because there are 
too few RCGs in these fields and very little differential ex- 
tinction. We assume that the slope of the reddening line, 
Ri, is the same within each region, but it may be differ- 
ent in different regions. For every region we had a collection 
of groups of bins, and the values of centroids of the RCGs: 
((I)rc, {V — I) rc) for each group. We fitted these data with 
Eq.lJ^J. All Ri measured by different methods are consistent 
with each other within their errors. In the following analy- 
sis we use the results by method [3] because they have the 
smallest scatter in the fitting. 

In Fig. [5] we show the distribution of centroids for re- 
gions (A), (B), (C) and (D) with the best fit lines. The best 
fit value Ri, the error a, and the standard deviation sdev 
for all regions, and their average values are given in table 
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Figure 2. Centroids (filled circle) of RCGs in the «J), (V - /)) 
CMD in the region (A), (B), (C) and (D) (from top to bottom) 
measured with method [3], as described in the text. The best fit 
line, parameters and the standard deviation (sdev) are also given 
in the figure. 



We also give the the mean Ri for all regions except (A), 
because region (A) has a significantly smaller Ri than other 
fields, and it has a very small error. In region (A) (Top panel 
of Fig. we did not use groups with (V — /}rc > 4 be- 
cause they are close to detection limit in V^-band, and their 
center in the CMD is shifted systematically to brighter I. 
We believe that groups with (V — I)rc < 4 are not affected 
by this effect because Ri does not change when this limit is 
reduced to smaller values. 

Our values of Ri are consistent with lUdalskil |2003) 
and are significantly smaller than the standard reddening 
law (~ 1.5). The measured Ri are shown as a function of 
the Galactic longitude (I) and latitude (b) in Fig.|3] We can 
see that Ri is slightly different from region to region, but 
there is no strong systematic dependence on the Galactic 
coordinates I or b. There might be a weak dependence in 
I, dRi/dl = -0.0102 ±0.0049. This trend gives the largest 
difference in Ri of 0.08 in regions (H) and (J). We think that 
this trend can be neglected in the following analysis because 
we take the scatter of 0.085 in Ri into account during the 
error estimation in reddening and extinctions. 

In the following analysis we use the mean value of the 
reddening slope for all regions i.e. Ri = 0.964 with sdev = 
0.085. The values of /q.rc, obtained by fitting Eq. (J^J with 
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Figure 3. The ratio of total to selective extinction Rj = 
A[/E(V — I) for each region as a function of Galactic coordi- 
nates. The values of Rj correspond to the middle of the error 
bars, and the average (Ri) is shown with a vertical line. 



Table 2. Extinction corrected / magnitude of RCGs Iq rc mea- 
sured with fixing Rj = 0.964 for each region. 



Region Iq,rc °7 



sdev 



1 o,rc 



°7 



0,RC 



A 
B 
C 
D 
E 
F 
G 
H 
1 
J 
K 

Note: /, 



14.625 
14.700 
14.808 
14.576 
14.551 
14.491 
14.453 
14.371 
14.808 
14.772 
14.529 



0.003 
0.005 
0.009 
0.004 
0.003 
0.003 
0.004 
0.016 
0.005 
0.006 
0.005 



0.034 
0.036 
0.042 
0.034 
0.031 
0.028 
0.033 
0.077 
0.041 
0.036 
0.026 



14.652 
14.728 
14.835 
14.603 
14.578 
14.518 
14.480 
14.398 
14.836 
14.799 
14.557 



0.024 
0.024 
0.025 
0.024 
0.024 
0.024 
0.024 
0.028 
0.024 
0.024 
0.024 



RC represent the value after the zero-point 



correction. 



the fixed slope Ri = 0.964 for each region, are given in 
table [5] and are plotted in Fig. [1] as a function of Galactic 
coordinates I and b. A clear evidence of the bar is apparent 
between —4 < I < 4 in the upper pan el of this figure, which 
is consistent with lStanek et all jl997t) . Note that regions (C) 
and (K) are shifted with respect to the pattern indicated by 
regions (G), (F), (E), (D), (A), (B); they are on the other 
side of the Galactic plane, i.e. in b > 0. 
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Table 1 . The ratio of total to selective extinction Rj for each region. 
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Figure 4. The constant -Zo.RC (cf. Eq. 0) is shown for each 
region as a function of Galactic coordinates adopting a fixed value 
Rl = 0.964. The values correspond to the middle of error bars. 
Note region (C) and (K) are fields with b > 0. 



4 EXTINCTION MAP 
4.1 Relative Reddening 

In this section we estimate the mean RCGs colours for each 
bin (V — I)rc, and we transform them to the relative red- 
dening, _B(V-/)rc = (V— -Orc — (V — I)o,kc, assuming the 
intrinsic colour of RCGs (V — I)o,rc = 1-0. Then we esti- 
mate the relative extinctions in V-band (Av,rc) and /-band 
(Ai,nc) by Eq. © and ©. 

In the following analysis we adopt the reddening line 
with the slope Ri = 0.964 and a constant Iq,rc as given in 
tabled Thanks to fixing Ri and Io.rc, we can get accurate 
(V — I)rc with a smaller number of RCGs, i.e. with high res- 
olution in space (small bin) and in reddening (small group). 
Furthermore we introduce a new indicator (V — I) a n which 
is the average colour of all stars in each bin, to represents 
the level of extinction in each bin to arrange bins into group. 
Because the number of all stars are much larger than RCGs, 
we to get higher spacial resolution. 

We divide each field into a new set of small " bins" with 
the size in the range 64 x 64 to 157 x 157 pixels, chosen so 



that each bin has ~ 200 stars. Next, we measure the average 
colour of all stars in each bin, {V — 7) a ii- We also measure 
the average colour of RCGs in each "bin", (V — I)rc, fol- 
lowing method [1] described in the previous section, with a 
radius of 0.5 mag, allowing the center of RCGs to lie only 
on the reddening line. The initial values of {V — 7)rc for 
this method are estimated by using the parallelogram along 
this reddening line. This may be fairly uncertain statisti- 
cally as there are only several RCGs per bin, but may suf- 
fer from less systematics than {V — I) a ii- We compare the 
two sets of colours in Fig. |S| (upper sequence) for the field 
BUL_SC22. There is good correlation between {V — I)rc 
and (V — /)au. Other fields have similar trends, but the 
slops are not always 1. This good correlation imply that 
(V — J) a n can be a good extinction indicator, as suggested 
bv lPopowski. Cook fc Beckerl (12003ft . 

We arrange bins in a given field in order of extinction 
by using {V — I) all. Then we combine these bins into groups 
from low extinction to high extinction until each group is 
filled by ~ 100 RCGs. Note that groups in §E]have ~ 1, 000 
RCGs. 

In Fig.|S|we show an example of a CMD for two groups 
of bins in BUL_SC22, one with low extinction (filled circles) 
and one with high extinction (open circles). We calculate 
the average colour of RCGs (V — I)rc,s for each group of 
bins using the method described in previous paragraph. As 
described above, thanks to fixing Ri and Io,rc, we can get 
accurate {V — I)nc,g, because the RCGs colour are very sim- 
ilar to the colours of red giants which are somewhat brighter 
or fainter than RCGs. 

There is a very good correlation between these colours 
and (V — I)rc obtained for single bins, as shown in Fig. El 
(lower sequence). In this figure the typical width of (V— I}rc 
of bins in each group can be seen by the gaps of (V — 7)Rc, g , 
about < 0.1 mag depending on the number density of stars. 
We adopt (V — I)rc,s as the mean colour of RCGs for bins 
in this group, because these are based on the large statistics 
of RCGs. 

Some bins were not included in any group for variety of 
reasons: (1) a bin has fewer than 20 stars, (2) a bin has a 
very bright blue star (mimics low extinction) , (3) a bin has a 
very bright red star (mimics high extinction). (4) because of 
a very high extinction a large fraction of Bulge red giants is 
below the detection limit, so these bins can not be grouped 
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Figure 6. The mean colour of RCGs (V — /)rc m each bin as 
a function of the mean colour of all stars in the bin ((V — /) a n, 
upper sequence) and the mean colour of RCGs in each group 
((V-/)Rc,g> lower sequence) in the field BULJ3C22. (V-I)nc,g 
is shifted by +0.5 mag for clarity. 



V-I 

Figure 5. The CMDs of two groups in BULJ3C22. Filled and 
open circles represent stars in a group with low and high ex- 
tinction, respectively. The elongated parallelogram represents the 
window used for the initial selection of RCGs. Two large circles 
indicate the RCG selection for the two groups. 



properly with {V-I) an (in fields BULJ3C5 and BUL_SC37). 
If there were more than 5 RCGs in a given bin, then we 
adopted the average colour of these RCGs (V — I)rc, which 
were estimated in the third paragraph of this section, as 
the mean colour of RCGs for this bin. Otherwise these bins 
are filled by the average value of the neighbours. {V — I)nc 
in each bin can not be used for the case (3), because stars 
reddened by bright red stars contaminate the RCG region 
in the CMD. 

Problematic bins are flagged with (1), "-2" (2), "- 
3" (3) and "-4" (4). If (V - J) RC is adopted in a bin then 
this bin is flagged with the above value plus "-10". 

We do not have confidence in measurements of {V— I)nc 
at (V — J)rc > 4, where the detection limit makes RCG 
centroid in the CMD systematically brighter and bluer. This 
implies that for (V — I)rc > 4 mag our maps give a lower 
limit to the extinction. We added "-20" to the value of the 
flag of such bins. Only BUL.SC5 and BULJ3C37 suffer from 
this effect. 



4.2 Zero point 

We assume that the average colour of RCG stars, cor- 
rected for interstellar reddening, is the same in every field, 
and we ignore a possible weak dependence on metalicity 
iPaczvnski fc Sta nek|l9 98j). The zero-po int of (V — I)o,nc is 
calibrated following lAlcock et al.llll998bl) . In Fig. [7] we show 
the offset of our relative extin ction Av.rc and V-ba nd abso- 
lute extinction for 20 RRab in lAlcock et all (|l998b), 
The average offset is (A v ,nc - ^4v,rr> = 0.055 ± 0.048 with 
the standard dev iation of sdev = 0.22 . The errors and sd ev 
are the same as in lAlcock et al.1 (Il998bl) with lStanekl jl996fi 's 
map. This means that reddening corrected (V— 7) colour and 
(J) magnitude of RCGs are given as 1 + 0.055/ Rvi = 1.028 
and the value of 7o,rc in table|2]plus 0.055Ri /Rvi ~ 0.027, 
respectively. The extinction corrected /-band magnitudes 
7q RC after zero-point correction are shown in table |2] 

Our extinction and reddening maps are calibrated by 
this offset. We show the Av extinction maps in Fig- [HI The 
parameters: the bin-size, the number of groups, the total 
number of all stars, the average number of all stars in each 
bin, the total number of RCGs, the average number of RCGs 
in each group, the average values and errors in E(V — /), 
Ay and Ai are given in Table |3] The values of (a) are the 
statistical errors in relative maps, and they do not include 
zero-point errors cte(v-i),o = 0.024, ua v ,o ~ 0.048 and 
o"A 7 ,o = 0.024. To check our maps we present in Fig. [5] a 
comparison of Av in the relatively large overlap region of 
BUL_SC30 and BUL_SC31. One can see a good correlation 
between them in this figure, and \ 2 /d.o.f. = 1.15 implies 
that our error estimate is reasonable. 
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Table 3. Basic parameters of maps, the bin-size, the number of groups, the total number of all stars, the mean number of all stars in 
each bin, the total number of RCGs, mean of RCGs in each group, the mean of values and errors in E(V — I), Ay and Aj. 
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Note: (a) do not include zero-point errors ce(v—i),0 = 0.024, OA v ,a = 0.048 and &Aj,o = 0.024. 



We show the histogram of E(V — I) for all our maps 
in Fig. 1101 The vertical line indicates the threshold value 
E(V — /)th = 3.0, which corresponds to {V — I)nc = 4 
mag. We have confidence in our maps below this threshold. 
However, small, very heavily obscured regions are above this 
reddening threshold, pushing some RCGs below the OGLE- 
II detection limit, and distorting the RCGs. This implies 
that for (V' — J)rc > 4 mag our maps give a lower limit to 
the extinction. Histograms of Av and Ai can be obtained 



by multiplying E(V — I) by Rvi and Ri respectively. The 
thresholds are: Ay.th = 5.8 and Ar,th = 2.9. 

The reddening corrected value of (V — /)o,rc may 
vary from one field to another due to weak depen- 
dence on metalicity. The RR Lyrae variables ar e good 
distance indicators (e.g., iNemec. Nemec fc Lutzl Il994t) 
whose period-luminosity relations ar e well established 
fljones. Carney. Storm fc Latharrl ll992b The RR Lyrae he 
in the instability strip, the range of their colour is small and 
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Figure 8. Ay Extinction map. Distant fields are in small boxes. Gray scale is shown at the right corner. 
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;A vrc~ A vrr > = - 055 (0.048), n=20, sdev= 0.21E 
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Figure 7. Difference of relative extinction .Ay nc from RCGs 
(see <I4. 1 1 and V-band extinction from 20 RRab in Alcock et al. 
(1998b), Ay r,Ri i.e., the same figure as Fig. 3 of Alcock et al. 
(1998b). 
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Figure 10. A histogram of E(V — I) for all our maps. Vertical 
line indicates the threshold value E(V — I)th = 3.0 mag. 




Figure 9. Comparison of Ay of overlap region in BUL_SC30 
and in BULJ3C31. 



it weakly depends on period, amplitude and/or metallicity 
jBono. fc StellingwerfM1994t lAlcock et alJll998afl . To check 
the relative zero-points of our reddening maps for each field 
we make use of the V — I colour of RR Lyrae Type ab 
(RRab), assuming that period - (V — I)o colour relation of 
RRab is the same for all OLGE-II fields 

We selecte d RRab in OGLE-II G alactic Bulge variable 
st ar catalogue llWozniak et al.ll200ll) by using the method 
of lAlard l dl996f) . We me asured periods P using PDM algo- 
rithm jstelhrigwerlll97^ ) and Fourier coefficients by fitting 
Fourier series with five harmonics. In Fig 1111 we show the 
the amplitude ratio R21 = A%jA\ versus phase differences 



02i = 4>2 ~ 20i for variables with P < 0.9 days, except 
0.4985 < P < 0.5001 days, which are affected by aliasing. 
We selected 1,961 RRab stars as a clear clump within the el- 
lipse in this figure. Of these 1,819 sta rs have (V) and (/) pho- 
tometry provided bv lUdalskil <l2003l) . We visually inspected 
all light curves and several stars with non RR Lyrae shape 
light curves were rejected. In Fig |12l we present colour - mag- 
nitude diagram for RRab stars so selected. The differential 
reddening is apparent. A similar reddening slope is found in 
ffl 

In the following analysis we rejected 39 stars with I < 
1.1(V — 7) + 13.1 because they are either nearby disk stars or 
they are blended with other bright stars. We did not reject 
several background RRab below the sequence as they made 
no difference to our analysis, for which we used 1,780 RRab 
stars. 

We assume that reddening corrected colour (V — 7)o,rr 
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V-I (mag) 

Figure 12. The CMD for RRab variables. Solid line with a slope 
of I/[V — I) = 1.0 is shown as a reference. RRab variables above 
the dashed line: / = 1.1(V — /) + 13.1 are rejected in our analysis 
because they are foreground disk stars or variables blended with 
other bright stars. 



is the same in every field, with a possible weak dependence 
on period P, or amplitude A, or both. We estimate the 
zero-points C% for each field i by fitting extinction corrected 
colour of RRab with a linear dependence on log P, A, and 
on both of them. We obtained the following relations: 

(V-I)o,rr = 0.309(±0.049)logP + C 4 ,sde« = 0.158(6) 

(V - J)o,rr = -0.096(±0.016)^ + d, sdev = 0.157, (7) 

(V-I)o,mt = 0.211(±0.056)logP- 0.063(±0.018)A 

+ d,sdev = 0.156. (8) 

All three relations provided almost the same values of 
d. Fitting with both P and A seems to be an over- 
parameterization. 

In Fig. 1131 we show extinction corrected colours of 
RCGs (V — J)o,rc (open circle) and RRab (V — 7)o,rr 
(filled circle and errors) as a function of the galactic I and b, 
where (V — 7)o,rc are measured by the method [3] in §[3]and 
(V — J)o,rr are the value given by Eq. © at logP = —0.3. 

(V — 7)o,rc is constant as assumed. The mean of 
(V — I)o,rc = 1.0283 ± 0.024, where we used a E (v-i),o = 
0.024 as the error in (V — 7)o,rc, is consistent with the 
colour of nearby RCGs (V — 7)o,RC,near = 1.01 ± 0.08 mag 
iPaczvriski fc Sta neklll99 8|). contrary to (V — J)q,rc bulge = 
1.11 ± 0.12 bv lPaczvhski et alJ Jl99Sl) with lStanekl lll996ft 's 
map. 

iy — 1)o,rr seems to be systematically redder by about 
~ 0.1 mag at large Such trend can not be seen in b. We 
plot (V — I)o,rr as a function |Z| in Fig. 1141 and we fit it 



with a straight line for \l\ < 6 (filled circle and solid line) 
and for \l\ > 6 (open circle and dashed line). One can see 
the RRab colour is constant at j/| < 6 and have a significant 
(3cr) dependence on \l\ at |/| > 6. The mean (V — J)o,rr is 
0.4670 ± 0.0047 with sdev = 0.0286 for \l\ < 6 and 0.5386 ± 

0. 0126 with sdev = 0.0773 for \l\ > 6. The significance of 
the redder colour in \l\ > 6 is estimated to be 5.3a level. 

This implies that (V — 7)o,rc, (V — I)o,rr or both of 
them, may vary with /. If we assume that (V — I)o,rr is 
constant, then the intrinsic colour of RCGs would be bluer, 

1. e., the zero-points of our reddening maps would be smaller 
by ~ 0.1 mag at \l\ — 10, compared to \l\ < 6. This could be 
explained if the outer RCGs have lower metalicity and are 
somewhat bluer than near the Galactic Center. The colour 
of RRab (V - 7)o,rr = 0.4 ~ 0.5 for \l\ < 6 is consis- 
tent with the colour of local RRab, while the two sets of 
co lours were claimed to be differen t bv lStanekl (^96) 's map 
bv lStutz. Popowski fc Gouldlll999l 

The scatter of (V — 7)o,rr values is large because public 
domain V^-band OGLE photometry is the average of ran- 
domly distributed small number of l/-band measurements, 
while the V amplitude is large. Furthermore the scatter is 
larger for \l\ > 6 because the number of RRab is small in 
these fields. Therefore, we do not correct for the offsets ap- 
parent in Fig. 1131 It will be possible to improve the accuracy 
of (V — 7)o,rr when individual V-band measurements be- 
come available. 



5 DISCUSSION AND CONCLUSION 

We confirmed the anomalous reddening, i.e. a small value of 
the ratio of total to selective extinction Ryi = 1.9 ~ 2.1 de - 
pending on the line of sight, as measured bv lUdalskiH2003l) . 
This implies that the distribution of dust grains may be 
tipped to smaller sizes in these regions, compared to aver- 
age in the Galaxy. A detailed analysis of these results is 
beyond the scope of the present study. 

By adopting the mean value of Ryi = 1.964, we have 
constructed reddening E(V — I) and extinctions, Ay and 
Ai maps in 48 OGLE-II GB fields, covering a range of 
— 11° < I < 11°, with the total area close to 11 square 
degrees. The reddening E(V — I) and extinctions, Ay and 
Ai are measured in the range 0.42 < E(V — I) < 3.5 mag, 
0.83 < Ay < 6.9 mag and 0.41 < Ai < 3.4. Note: above 
the threshold values E(V — I)th = 3.0, Ay tt h = 5.8 and 
A/. tn = 2.9, which correspond to (V — I)nc = 4 mag, our 
maps give the lower limit to reddening and extinction. Spa- 
tial resolutions of maps are 26.7" ~ 106.8" depending on 
the stellar density of each field. 

The absolute zero-point is calibrated us ing 20 R Rab 
variable s in Baade's Window, following lAlcock et alJ 
J19981J) . Relative zero-points of our maps are verified with 
1,780 RRab variables found in our fields. We found that 
these zero-points may be lower by ~ 0.1 mag at larger Galac- 
tic longitudes \l\ > 6. Note that our extinction map is OK 
in terms of the relative extinction within each field. The rel- 
ative zero points are also OK between fields for |/| < 6. We 
did not make any correction for this effect in the paper. 

We used the mean value of Ryi = 1.964, but differ- 
ent fields have a range of values, with a standard deviation 
sdev ~ 0.085 (see Figl^J. We estimated the errors of our 
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Figure 13. Extinction corrected colours of RCGs (V — /)o,R,c 
(open circle) and RRab (V — /)o,rr (filled circle and errors) as 
a function of the galactic I and b. (V — I)o rr are the value 
given by Eq. JSJ at logP = —0.3. Solid line represent a mean of 
(V - I)o,B.C = 1-0283 ± 0.024, sdev = 0.0034, where the error in 
(V - I)o,nc is dominated by <t E (V-I),q 



: 0.024. 



maps by taking this range into account. The errors approach 
0.17 mag at the highest extinct i on (A y = 6 ). 

As noted bv lUdalski et alJ (120021) and lUdalskil J2003ft . 
the /-band filter used by OGLE-II has the red wing some- 
what wider than standard. This may lead to systematic de- 
viations from the standard values for very red stars, giv- 
ing brighter /-band magnitudes (redder V — I colour) for 
the OGLE-II filter for very red stars (V - I > 2), while 
the error is negligible in the range where the OGLE-II 
data were calibrated by standards, i.e. for (V — I < 2). 
This effect may reach ~ 0.2 mag for the very red RCGs 
(V — I > 4 mag). This effect makes the slope of re ddening 
line in the CMD somewhat shallower, as shown bv lUdalska 
(2003). This effect is not sufficie nt to explain th e anoma- 
lous extinction towards the GC. lUdalskj ([2003) analyzed 
the differe nces between O GLE-II and standard /-band fil- 
ters using iKuruczl |l992) 's model atmosphere of a typical 
R CGs, reddened with the standard int erstellar extinctio n 
of lCardelli. Clavtonfc Mathisl Jl989l) and lFitzpatrickl J1999D . 
He found possible errors to be at the level ±0.1 in Rvi, de- 
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(V-I) 0EE = 0.4662(0.0077) + 0.0004(0.0029) 1, 
X 2 /(35_ 2 )=l 00. sdev=0.029 
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Figure 14. Extinction corrected colours of RRab (V — /)o,rr as 
a function of the absolute galactic \l\. (V — I)q rr are the value 
given by Eq. JHJ at logP = -0.3. Filled and open circles represent 
plots in |Z| < 6 and in / > 6, and solid and dashed lines indicate 
the best fits for them, respectively. One can see the RRab colour 
is constant at |2[ < 6 and significantly red at |/| > 6. 



pending on the model of extinction and the colour range of 
RCGs in each field. This leads to ~ 0.2 mag differences in 
E(V-I), A v and Ai for the very red RCGs, with V-I > 4 
mag, corresponding to E(V~ I) ~ 3, Ay ~ 5.9 and Aj — 2.9 
mag in our extinction maps. These differences are the up- 
per limits because the atmospheric absorption in the range 
900-990 nm m akes the OGLE-II / -band filter closer to the 
standard one llUdalski et alJ 12002). To make a more accu- 
rate estimate of the required correction it is necessary to 
have very red standard stars, with (V — /) > 2. 

The extinction corrected /-band magnitude of RCGs, 
Iq rc = 14.6 in Baade's window is fainter than expected. 
Adopting distance mod ulus to the GC of 14.52 ± 0.1 mag 
jEisenhauer et alJl2003T) . and assuming that the population 
of RCG stars in the Galactic bulge is the same as local, 
i.e. that the absol ute magnitude is Mi RC = —0.26 ± 0.03 
iAlves et alJ |20 02[) and the average co lour is (V — I)o = 
1.01 ± 0.08 iPaczvriski fc Staneklfl99ci) . then the expected 
magnitude of RCG in Baade's Window should be Io,rc = 
14.3. In Fig. EH1 we show CMD and RCG centroids of 
BUL_SC1 (at Baade's window) before (left panel) and af- 
ter (right panel) the extinction and reddening correction, 
with the expected RCG centroid. Note, that the difference 
between the average distance modulus in Baade's Window 
and the Gal actic center is only 0.02 m ag, i.e. it is of no 
consequence l|Paczvriski fc St anck 1998J). 

We do not know what is the solution of this problem. 
It may be that the population ef fects are larg e Mi n or , — 
0,4 ~ -0.03, as claim ed by ISalaris et aJT i2003l) and 
iPercival fc Salaris I ll2003h . it may be that the distance to the 
GC is larger, or it is possible that the reddening is more com- 
plicated. We assumed that the extinction to reddening ratio 
is constant all the way to zero extinction. OGLE photometry 
is well calibrated with standards for (V — /) < 2. The slope 
of the reddening line in Baade's Window (Region D) is well 
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V-I (mag) (V-I) (mag) 

Figure 15. Colour Magnitude Diagram of BUL_SC1 before (left panel) and after (right panel) the extinction and reddening correction. 
Small and large filled circles indicate the same plots as Fig. [5] and RCGs centroid after the extinction and reddening correction. Open 
circle represents the RCGs centroid expected from the assumption that the population of RCGs in GB is the same as local. The solid 
line indicate the reddening line given by Eq. [3]with the mean slope Rj = 0.964 and a constant Iq rq = 14.6 for this field. 



measured for RCG in the colour range 1.6 < (V — I) < 2.2, 
but we have no direct information about the reddening 
line for 1.0 < (V — I) < 1.6, i.e. for the reddening range 
< Ev-i < 0.6. If we make an ad hoc assumption that the 
RCG population in Baade's Window is the same as local, 
and the distance is 8 kpc, then we may obtain Iq ro = 14.3 
and (V-I) ,rc = 1.0 adopting Rvi ~ 2.8 (i.e. 7?/ ~ 1.8)for 
the unobserved range < Ev-i < 0.6. We do not know if 
this is plausible or not. 

The situation will improve somewhat once detailed V- 
band OGLE photometry becomes available for RR Lyrae 
stars, and this will make it possible to obtain an indepen- 
dent estimate of the reddening in Baade's Window. Prelim- 
inary analysis of the average photometry of RR Lyrae stars 
seems to indicate that the puzzling brightness of red clump 
giants may be due to population effects, i.e. their average 
absolute magnitude in the Bulge is somewhat different than 
it is near the Sun. A much improved analysis will be done 
when OGLE-II V^-band measurements become available in 
the next several months. At this time it should be OK to use 
our maps of differential reddening, but we consider our cal- 
ibration of zero point to be uncertain. This implies that at 
this time our reddening maps are not adequate for a quan- 
titative study of Galactic bar structure. 

In this work we leave our extinction maps in the OGLE- 
II /-band because the reddening map is urgently needed 
for various applications, and is already used in some works 



JSumi et all2003rilWrav. Ever fc Paczvhskil; 2003L A reader 
must take care of possible small error mentioned above while 
using our maps for standard /-band photometry. 

The extinction map is available in elec- 
tronic format via anonymous ftp from the server 
ftp:/ /ftp. astrouw. edu.pl/ ogle/ ogle2 / extinction/ and 
ftp://bulge.princeton.edu/ogle/ogle2/extinction/ 

These extinction maps of OGLE-II GB fields facilitate 
the study of th e Galactic structur e with OGLE proper mo- 
tion catalogue (Sumi et al. 2003b) and microlensing optical 
depth, and a study of variable stars, but the reader should 
be aware of the zero point of the extinction may not be accu- 
rate. We intend to improve the quality of the zero points in 
all fields as soon as individual OGLE V-band measurements 
of the RR Lyrae stars become available. 
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